Blackbody radiation shift in the Sr optical atomic clock 
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We evaluated the static and dynamic polarizabilities of the 5s 2 1 So and 5s5p 3 Pq states of Sr using 
the high-precision relativistic configuration interaction + all-order method. Our calculation explains 
the discrepancy between the recent experimental 5s 2 1 So — 5s5p 3 Pq dc Stark shift measurement 
Aq = 247.374(7) [Middelmann et. al, larXiv:1208. 28481 (2012)] and the earlier theoretical result of 
261(4) a.u. [Porsev and Derevianko, Phys. Rev. A 74, 020502R (2006)]. Our present value of 
247.5 a.u. is in excellent agreement with the experimental result. We also evaluated the dynamic 
correction to the BBR shift with 1% uncertainty; -0.1492(16) Hz. The dynamic correction to the 
BBR shift is unusually large in the case of Sr (7%) and it enters significantly into the uncertainty 
budget of the Sr optical lattice clock. We suggest future experiments that could further reduce the 
present uncertainties. 
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I. INTRODUCTION 

Optical lattice clocks have shown tremendous progress 
in recent years An optical frequency standard based 
on the 5s 2 1 So — 5s5p 3 Pg transition of ultracold 87 Sr 
atoms confined in a one-dimensional optical lattice is pur- 
sued by a number of groups 043 • Its systematic uncer- 
tainty has been demonstrated at the 10~ 16 fractional fre- 
quency level and an order-of-magnitude improvement is 
expected to be achieved soon [J [j] . A three-dimensional 
optical lattice clock with bosonic 88 Sr was demonstrated 
for the first time in 0. 

The measured clock transition frequencies must be cor- 
rected in practice for the effect of the ambient blackbody 
radiation (BBR) shift, which is quite difficult to measure 
directly. The BBR shift can only be suppressed by cool- 
ing the clock. At room temperature, the differential BBR 
shift of the two levels of a clock transition turns out to 
make one of the largest irreducible contributions to the 
uncertainty budget of optical atomic clocks. The Sr clock 
transition has the largest BBR shift of all optical fre- 
quency standards that are currently under development 
(see Ref. @ for a recent review). The fractional BBR 
shift A^bbr/^o in Sr is more than a factor of 1000 larger 
than the fractional BBR shift in the Al + ion clock [lfj. 
The BBR shift of an optical clock can generally be ap- 
proximated by the dc Stark shift of the clock transition 
to about 1-2% precision, because optical frequencies are 
100 times greater than characteristic BBR frequencies. 
However, Sr represents an exception, where the so-called 
dynamic correction [lit] , that needs to be determined sep- 
arately from the dc Stark shift, is 7%. Recently, the dc 
Stark shift in Sr has been measured with 0.003% pre- 



cision [131 ]. and the dynamic correction was evaluated 
based on a set of El transition rates and the Stark shift 
measurement. The measured value differed substantially 
(by almost Act) from the previous theoretical determina- 
tion [ll|. 

In this work, we evaluate the static and dynamic po- 
larizabilities of the 5s 2 1 So and 5s5p 3 Pq states of Sr us- 
ing the high-precision relativistic CI+all-order method. 
Our calculation explains the above-mentioned discrep- 
ancy between the experimental 5s So — 5s5p 3 Pq dc 
Stark shift measurement Aa = 247.374(7) a.u. [l3| ] and 
the earlier theoretical result of 261(4) a.u. PJJ. We found 
that the El matrix elements for the transitions that give 
dominant contributions to the 3 Pq polarizability, in par- 
ticular the 5s4d 3 Di — 5s5p 3 Pg, are rather sensitive to 
the higher-order corrections to the wave functions and 
other corrections to the matrix elements beyond the ran- 
dom phase approximation. A correction of only 2.4% to 
the dominant 3 Di— 3 Pq matrix element leads to 5% dif- 
ference in the final value of the 3 Pg — 1 So Stark shift. 
In this work, we included the higher-order corrections in 
an ah initio way using the CI+all-order approach, and 
also calculated several other corrections omitted in [ill ]. 
Our value for the dc Stark shift of the clock transition, 
247.5 a.u., is in excellent agreement with the experimen- 
tal result 247.374(7) a.u. [lj. 

We have combined our theoretical calculations with the 
experimental measurements of the Stark shift [131 ] and 
magic wavelength Q of the 5s 2 1 So — 5s5p 3 Pq transition 
to infer recommended values of the several electric-dipole 
matrix elements that give the dominant contributions to 
the 3 Pq polarizability. We used these values to evaluate 
the dynamic correction to the BBR shift of the """So — 3 Pq 
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TABLE I: Comparison of experimental P3] and theoretical energy levels of Sr in cm -1 . Two-electron binding energies are 
given in the first row, energies in other rows are given relative to the ground state. Results of the CI+MBPT and Cl+all- 
order calculations are given in columns labeled "CI+MBPT" and "CI+A11 (A)". The Cl+all-order values with the ground 
state two-electron binding energy shifted by 200 cm -1 are given in column labeled "CI+A11 (B)". Corresponding relative 
differences of these three calculations with experiment are given in the three corresponding columns labeled "Diff." in %. The 
5s4d 3 Di — 5s5p 3 P§ transition energy is given in the last row. 
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transition to be -0.1492(16) Hz. 

We determined that the 5s4d 3 Di — 5s5p 3 Pg tran- 
sition contributed 98.2% to the dynamic correction for 
the 3 Pq level. Our calculation enables us to propose 
an approach for further reduction of the uncertainty in 
the BBR shift. In particular, there is a correlation in 
the uncertainty of the BBR shift and the lifetime of the 
hsAd 3 Di state, if branching ratios are known to sufficient 
accuracy. At present, experimental measurements of the 
5s4d 3 D term-averaged lifetime have an uncertainty of 
about 7% [3, HEj . We note that the experiment [15j . 
which was performed at JILA some 20 years ago, has rel- 
evance in the determination of the uncertainty budget of 
one of the world's most accurate clocks now being devel- 
oped at the same institution - a development probably 
not envisaged at the time. 

A new determination of this (or 3 Di) lifetime with 
0.5% uncertainty would provide a value of the Sr clock 
BBR shift that is accurate to about 0.5%, which would 
be a factor of 2 improvement in the uncertainty that 



we state here. This result is determined by the rele- 
vant branching ratios needed for the extraction of the 
5s4d 3 Di — 5s5p 3 Pg matrix elements from the lifetime 
measurement. We have determined these branching ra- 
tios with an uncertainty of 0.2%. A further reduction in 
the uncertainty of the Sr clock BBR shift could be ef- 
fected by an improved measurement of these branching 
ratios. The lifetime of the corresponding 6s5d 3 T>i state 
in Yb has been recently measured in Ref. [16J. 



II. METHOD AND ENERGY LEVELS 

Calculation of Sr properties requires an accurate all- 
order treatment of electron correlations. This can be 
accomplished within the framework of the Cl+all-order 
method that combines co nfig uration interaction and 
coupled-cluster approaches [HHHI ■ To evaluate un- 



certainties of the final results, we also carry out CI 
and Cl+many-body perturbation theory (MBPT) 



22] 
23] 
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TABLE II: The CI+MBPT and Cl+all-order results and further corrections to the El matrix elements for transitions that 
give dominant contributions to the polarizabilities of the 5s 2 1 So and 5s5p 3 Pq states. The CI+MBPT and Cl+all-order 
results including RPA corrections are given in columns labeled "MBPT+RPA" and "All+RPA", respectively. The relative 
differences between the CI+all-order+RPA and CI+MBPT+RPA results are given in column "Higher orders" in %. The 
other contributions include the core-Brueckner (<r), two-particle (2P), structural radiation (SR), and normalization (Norm) 
corrections. Total relative size of corrections beyond CI+all-order+RPA is given in column "Corr." in %. The recommended 
values for the 5s 2 1 So — 5s5p 1 P° matrix element was obtained from the 1 P° lifetime measurement [13], and the recommended 
values for all other transitions are from the present work (see Section |V)) . 
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SR 
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Corr.(%) 
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-1.23 
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-0.016 
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-1.69 
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5s5p 3 Pg 
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calculations. These methods have been described in a 
number of papers [l^, [H, [H, HI] and we provide only a 
brief outline of these approaches. 

We start with a solution of the Dirac-Fock (DF) equa- 
tions 

H ip c = e c i> c , 

where H is the relativistic DF Hamiltonian [l8], [23| and 
ip c and e c are single-electron wave functions and energies. 
The calculations are carried out in the V N ~ 2 potential. 
The wave functions and the low-lying energy levels are 
determined by solving the multiparticle relativistic equa- 
tion for two valence electrons [221 ] . 

H eS (E n )<Z> n = E n $ n . 

The effective Hamiltonian is defined as 

H eG (E) = H FC + £(£), 

where Hfq is the Hamiltonian in the frozen-core approx- 
imation. The energy-dependent operator £(£") which 
takes into account virtual core excitations is constructed 
using second-order perturbation theory in the CI+MBPT 
method [23| and using a linearized coupled-cluster single- 
double method in the Cl+all-order approach 18]. It is 
zero in a pure CI calculation. We refer the reader to 
Refs. [HI, [23| for detailed description of the construction 
of the effective Hamiltonian. 

Unless stated otherwise, we use atomic units (a.u.) for 
all matrix elements and polarizabilities throughout this 
paper: the numerical values of the elementary charge, 
e|, the reduced Planck constant, Ti — h/2n, and the 
electron mass, m e , are set equal to 1. The atomic 
unit for polarizability can be converted to SI units via 
a/h [Hz/(V/m) 2 ]=2.48832xl0" 8 a (a.u.), where the con- 
version coefficient is 4ireociQ/h and the Planck constant 
h is factored out in order to provide direct conversion 
into frequency units; Oo is the Bohr radius and eo is the 
electric constant. 

As a first test of the accuracy of our calculations, we 
compare our theoretical energies with experiment for a 



number of the even and odd parity states. Comparison of 
the energy levels (in cm" 1 ) obtained in the CI+MBPT, 
and Cl+all-order approximations with experimental val- 
ues [IH is given in Table HI Ground state two-electron 
binding energies are given in the first row of Table HI en- 
ergies in other rows are measured from the ground state. 
The relative differences of the CI+MBPT and Cl+all- 
order calculations with experiment (in %) are given in 
columns labeled "DifF". Since the Cl+all-order values 
are systematically higher than the experimental values, 
a large fraction of the difference from experiment can be 
attributed to the difference in the value of the ground 
state two-electron binding energy. We find that shifting 
the Cl+all-order value of the ground state two-electron 
binding energy by only 200 cm -1 (see results in column 
Cl+all-order (B)) brings the results into excellent agree- 
ment with experiment for most of the states. We give the 
hsAd 3 Di — bs5p 3 Pq transition energy in the last row of 
Table [I] This transition is particulary important to the 
subject of this work, since it contributes 61% to the static 
polarizability and 98% to the dynamic correction to the 
BBR shift of the 5s5p 3 Pq state. In fact, the acciden- 
tally small value of this transition energy is the source 
of the anomalously large (7%) dynamic correction to the 
BBR shift of the ^Sq -+ 3 P ° transition in Sr. We see con- 
siderable improvement of the accuracy in this transition 
energy from the CI+MBPT to Cl+all-order approxima- 
tion, by a factor of 6. The CI+MBPT and Cl+all-order 
values differs from the experiment by 11% and 1.7%, re- 
spectively. 

III. AB INITIO CALCULATION OF 
ELECTRIC-DIPOLE MATRIX ELEMENTS 

The reduced electric-dipole matrix elements are ob- 
tained with the Cl+all-order wave functions and effec- 
tive electric-dipole operator D c g in the random-phase 
approximation (RPA). The effective operator accounts 
for the core-valence correlations in analogy with the ef- 
fective Hamiltonian [24], HBJ • We include additional cor- 
rections beyond RPA in the calculation of the El matrix 
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TABLE III: Contributions to the 5s 2 1 So and 5s5p 3 Pq static polarizabilities of Sr in a.u. The dominant contributions to the 
valence polarizabilities are listed separately with the corresponding absolute values of electric-dipole reduced matrix elements 
given in columns labeled D. The theoretical and experimental 12] transition energies are given in columns AE t h and A_E e x P t- 
The remaining contributions to valence polarizability are given in rows Other. The contributions from the core and a vc terms 
are listed together in rows Core + Vc. The dominant contributions to ao, listed in columns «o[A] and «o(B), are calculated 
with CI + all-order +RPA (no other corrections) matrix elements and theoretical [A] and experimental [B] energies [12| ]. 
respectively. The dominant contributions to ao listed in column ao [C] are calculated with experimental energies and our final 
ab initio matrix elements. The dc Stark shift for the 5s5p 3 Pq — 5s 2 1 So transition is listed in the last rows of the table. 



State 


Contribution 
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^ r a l 
ao [A] 


a [BJ 
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~ rr 1 ! 
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5.208 
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0.25 


0.25 




0.25 
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0.95 


0.95 




0.95 




Other 








4.60 


4.60 




4.60 




Core + Vc 








5.29 


5.29 




5.29 




total 








19/ .8 


1 HQ CI 




ly4.4 




Recomm.' c ' 














19( .14(20) 


5s5p 3 Pg 


5s5p 3 P - 5s4d 3 Di 


3777 


3842 


2.712 


285.0 
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14721 


1.970 


38.7 


38.6 


1.940 


37.4 




5s5p 3 Pg - 5s5d 3 Di 


20660 


20689 


2.460 


42.9 


42.8 


2.432 


41.8 




5s5p 3 Pg - 5p 2 3 Pi 


21208 


21083 


2.619 


47.3 


47.6 


2.620 


47.6 




5s5p 3 Pg - 5s7s 3 Si 


23056 


23107 


0.516 


1.69 


1.69 




1.69 
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1.161 


7.8 


7.8 




7.8 
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29.1 


29.1 




29.1 
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5.55 




5.55 




Total 
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247.374(7) 



a CI+all-order+RPA values (no other corrections). b CI+all-order+RPA + other corrections. c Obtained using experimental 
5s5p X P? lifetime from d Obtained using recommended value for the 1 So polarizability and the experimental value of the 

Stark shift [l]|. 



elements in comparison with |ll|, |26| . These contribu- 
tions include the core-Brueckner (u), two-particle (2P) 
corrections, structural radiation (SR), and normalization 
(Norm) corrections [U, [25[ . While we find some cance- 
lation between the various corrections, these cannot be 
omitted at the 1% level of accuracy. Partial cancela- 
tion of the structural radiation and normalization cor- 
rections was discussed in Ref. [27|. Detailed analysis of 
the structure radiation correction was carried out in the 
same work (27| . 

The results for several transitions that give dominant 
contributions to the 1 So— 3 Pg dc Stark shift are sum- 
marized in Table |TTJ The percentage differences between 
the CI+all-order+RPA and CI+MBPT+RPA calcula- 
tions are given in the column labeled "Higher orders" . 
We note that it is positive for some transitions and neg- 
ative for other transitions. Our final ab initio values are 
given in column labeled "Final" . We find that total rela- 
tive size of corrections beyond CI+all+RPA given in col- 
umn labeled "Corr" is small, 0.04-1.7%, but significant. 
We estimate the uncertainties in the ab initio values of 
the matrix elements to be 1% based on the comparison 



of the CI+MBPT+RPA and CI+all-order+RPA values 
and combined size of other corrections. 

We also provide the recommended values for these 
transitions. The recommended value for the 5s 2 1 So — 
5s5p X P° matrix element was obtained in [HlUfil from the 
1 P 1 ° lifetime measurement from photoassociation spectra 
[171 ] , the recommended values for all other transitions are 
obtained in the present work in Section [V] 



IV. POLARIZABILITIES 

We evaluated the static and dynamic polarizabilities 
of the 5s 2 1 So and 5s5p 3 Pg states of Sr using the high- 
precision relativistic Cl+all-order method. The scalar 
polarizability ao(ui) is separated into a valence polariz- 
ability Oq(u>), ionic core polarizability a c , and a small 
term a vc that modifies ionic core polarizability due to 
the presence of two valence electrons. The valence part 
of the polarizability is determined by solving the inhomo- 
geneous equation in valence space, which is approximated 



TABLE IV: Breakdown of the contributions to the 5s5p 3 Pq 
static polarizability ao(uj — 0) and dynamic polarizability 
«o(w) at the 813.4 nm magic wavelength. The dominant con- 
tributions to the valence polarizabilities are obtained with 
experimental energies and recommended values of the ma- 
trix elements. The electric-dipole reduced matrix elements 
are given in column labeled "j) rccom " The experimental [12] 
transition energies are given in column labeled "AU cxpt " . The 
remaining contributions to valence polarizability are given in 
row labeled "Other". The contributions from the core and 
a vc terms are listed together in row labeled "Core + Vc". 
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FIG. 1: (Color online) The frequency-dependent polarizabil- 
ities of the Sr 5s 2 1 So and 5s5p 3 Pg states near 813.4 nm 
magic wavelength. The frequency-dependent polarizability of 
5s5p 3 Pg shifted by ± 1% are shown to illustrate the sensi- 
tivity of the magic wavelength to the 5s5p 3 Pg polarizability. 
The magic wavelength is marked with arrow. 



as m 



(E v -H eS )\m{v,M')) =D eS \V (v, J,M)) (1) 



for the state v with total angular momentum J and pro- 
jection M. The wave function ^/(v,M') is composed 
of parts that have angular momenta of J' = J, J ± 1 
that allows us to determine the scalar and tensor polar- 
izability of the state \v,J,M) (28|. The effective dipole 
operator D c g includes RPA corrections. The core and 
a vc terms are evaluated in the random-phase approxi- 
mation. Their uncertainty is determined by comparing 
the DF and RPA values. The small a vc term is calcu- 
lated by adding a vc contributions from the individual 
electrons, i.e. a vc (5s 2 ) = 2a vc (5s), and a vc (5s5p) = 
a vc (5s) +a vc (5p). The frequency dependence of the core 
and a vc terms is negligible, and we use their static values 
in all calculations. 

While we do not use the sum-over-states approach in 
the calculation of the polarizabilities, it is important to 
establish the dominant contributions to the final values. 
We combine the electric-dipole matrix elements and en- 
ergies according to the sum-over-states formula for the 
valence polarizability [29j : 



3(2 J + 1) 



E 



(E n ~ E v )\(v\\D\\n)\< 

(E n - E v ) 2 - W 2 



(2) 



to calculate the contributions of specific transitions. 
Here, J is the total angular momentum of the state v, 
D is the electric-dipole operator, Ei is the energy of the 
state z, and frequency uj is zero in the static polarizability 
calculations. 

We have carried out several calculations of the dom- 
inant contributions to the polarizabilities using differ- 
ent sets of the energies and El matrix elements in or- 



der to understand the difference of the theoretical pre- 
dictions for the Stark shift 5s 2 1 S U — 5s5p 3 P ° Aa = 
261(4) a.u. and recent experimental measurement Aa = 
247.374(7) a.u. as well as to provide a recommended 
value for the 5s5p 3 Pg — 5s4d 3 Di matrix element. The 
results are summarized in Table Hill Other theoretical cal- 
culations of Sr polarizabilities were recently compiled in 
review [29] . The ground-state polarizability of Sr was cal- 
culated using relativistic coupled-cluster (RCC) method 
in 30]. Their value 199.7(7.3) a.u. is in agreement with 
our calculations. 



In Table IIIII the absolute values of the correspond- 
ing reduced electric-dipole matrix elements are listed in 
columns labeled "D" in a.u.. The theoretical and ex- 
perimental [12j transition energies are given in columns 
Ai? t h and A_E oxpt . The remaining valence contributions 
are given in rows labeled "Other". The contributions 
from the core and a vc terms are listed together in row 
labeled "Core +Vc". The dominant contributions to ao 
listed in columns ao[A] and ao[B] are calculated with 
CI + all-order +RPA (no other corrections) matrix ele- 
ments and theoretical [A] and experimental [B] energies 
[13], respectively. 

Our ao[A] result agrees with the earlier calculation of 
[ill ] which was carried out using CI+MBPT approach 
with energy fitting that approximated missing higher- 
order corrections to the wave functions. We note that 
this may be fortuitous since the calculation of [ll| was 
carried out in V N potential, while we are using V N ~ 2 
potential since the present version of the CI+all-order 
method is formulated for V N ~ 2 potential. The El ma- 
trix elements in [TJ, [2(| included RPA but omitted all 
other corrections calculated in the present work. We find 
that replacing the theoretical energies with experimental 
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TABLE V: Dynamic corrections to the BBR shift of the 5s 2 1 So - 5s5p 3 Pg clock transition in Sr at T = 300 K (in Hz). 
Quantities rji, AEg yn , and A^3p„_ ls are defined in text. 







??2 


V3 


V 


a (uj = 0) 


AE^ yn /h 


Total (5s 2 ^o) 


0.00163 


0.00001 





0.00164 


197.1 


-0.0028 


5s5p 3 Pq - 5s4d 3 Di 


0.03394 


0.00414 


0.00088 


0.03896 






5s5p 3 Po - 5s6s 3 Si 


0.00032 








0.00032 






5s5p 3 Pg - 5s5d 3 Di 


0.00018 








0.00018 






5s5p 3 Pq - 5p 2 3 Pi 


0.00019 








0.00020 






5s5p 3 Pg - 5s7s 3 Si 


0.00001 








0.00001 






5s5p 3 Pg - 5s6d 3 Di 


0.00002 








0.00002 






Total (5s5p 3 Pg) 


0.03467 


0.00415 


0.00088 


0.03970 


444.6 


-0.1520 


Final A^_ 1Sq 
Ref. [13] 












-0.1492(16) 
-0.1477(23) 



values reduces the Stark shift by 2.3%. We note that in 
the case of Sr all of the states contributing to the po- 
larizabilities are included in our computational basis and 
this procedure is not expected to cause problems with 
basis set completeness, as in the case of Yb 31]. The 
dominant contributions to cto listed in column ao[C] are 
calculated with experimental energies and final ab initio 
matrix elements. Inclusion of the small corrections fur- 
ther reduces the value of the Stark shift by 3.1%, and our 
resulting value obtained with our final ab initio matrix 
elements is in excellent agreement with experiment (T3 |. 



V. DETERMINATION OF RECOMMENDED 
VALUES OF ELECTRIC-DIPOLE MATRIX 
ELEMENTS 

To further improve our values of the other matrix ele- 
ments, we use the measurement of the magic wavelength 
for the 3 Pq — 1 So clock transition to determine recom- 
mended values of the 5s5p 3 Pq — 5s6s 3 Si, 5s5p 3 Pg — 
5s5d 3 Di, and 5s5p 3 Pg — 5p 2 3 P\ matrix elements. The 
polarizability at the 813.4 nm magic wavelength is 
286.0 a.u. which is essentially fixed by the value of the 
5s5p X P? lifetime H<|. The contributions to the 3 Pg 
polarizability at the magic wavelength are listed in Ta- 
bic [IVJ Since the contribution of the 5s5p 3 Pq — 5s6s 3 Si 
transition is dominant, the magic wavelength limits the 
value of this matrix element within about 0.5%. Since we 
appear to systematically overestimate the correction to 
matrix elements beyond CI+all+RPA approximation, we 
adjust the values of two other 5s5p 3 Pfj — 5s5d 3 Di and 
5sbp 3 Pq — bp 2 Pi matrix elements in a similar way as 
the 5s5p 3 Pq — 5s6s 3 Si one. We plot the dynamic polar- 
izabilities of the 1 So and 3 Pq states in the vicinity of the 
magic wavelength on Fig. Q] to illustrate that the cross- 
ing point is extremely sensitive to the matrix element 
values. The 0.5% change in the values of the matrix el- 
ements (corresponding to 1% change in the value of the 



3 -Pq polarizability) shifts the crossing point by more than 
4 nm. 

After we determined the values of these three ma- 
trix elements, we used them to obtain the value of the 
5s5p 3 Pq — 5s4d 3 Di matrix element from the experi- 
mental value of the Stark shift. We determine its un- 
certainty from the uncertainty of all the other contri- 
butions to the 3 Pq polarizability value (listed in the 
last column of Table II VI) . Since our theoretical values 
may experience a systematic shift in one direction, we 
(somewhat conservatively) simply add all of the uncer- 
tainties, totaling to 3.3 a.u, instead of adding them in 
quadrature. Assigning this value to be the uncertainty 
in the dominant 3 D\ contribution of 272.7 a.u., we esti- 
mate the uncertainty in the recommended value of the 
5s5p 3 Pg — 5sid 3 Di matrix element to be 0.5%. Since 
the contributions to both static and dynamic polarizabil- 
ities from 5s5p 3 Pg - 5s7s 3 Si and 5s5p 3 P£ - 5s6d 3 Di 
transitions are small, we use ab initio CI+all+RPA val- 
ues and assign them 1.5% uncertainty. Combined with 
experimental energies and other small contributions, the 
set of recommended matrix elements reproduces recom- 
mended values for both the 5s5p 3 Pq static and dynamic 
polarizability at 813.4 nm magic wavelength. 



VI. BLACKBODY RADIATION SHIFT 

The leading contribution to the multipolar black body 
radiation (BBR) shift of the energy level g can be ex- 
pressed in terms of the electric dipole transition matrix 
elements [HJ 

^E g = -^^-Y.M\D\\n)\ 2 F x {y n ). (3) 

Here ks is the Boltzmann constant, y n = (E n — 
Eg)/ (ksT), and F\ (y ) is the function introduced by Far- 
ley and Wing in [321 ] . Its asymptotic expansion is given 
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TABLE VI: Experimental transition energies (in cm -1 ), theoretical line strengths (in a.u.), transition rates (in s _1 ), and 
branching ratios for transitions contributing to the 5s4d 3 Dj lifetimes. The CI+RPA, CI+MBPT+RPA, and CI+all- 
order+RPA results are listed in columns labeled "CI", "MBPT", and "All", respectively. The recommended values of the 
(5s5p 3 P? i2 ||.D||5s4d 3 Di, 2 }| and |(5s5p 3 Pi||D||5s4d 3 D 3 )| matrix elements are obtained using the recommended matrix ele- 
ment for the 5s4d 3 Di — > 5s5p 3 Pg transition (i.e. scaled by 0.9862). Experimental energies are used in all cases. Numbers in 
square brackets represent powers of 10. 



Transition 






Line strengths S 




Transition rates A a b 


Branching ratios 






CI 


MBPT 


All 


Recomm. 


MBPT 


All 


Recomm. 


CI 


MBPT 


All 


3 D a -)■ A P° 


3842 


9.503 


7.189 


7.357 


7.156 


2.753(5] 


2.817(5] 


2.740(5] 


0.5949 


0.5954 


0.5953 


3 Di 3 P? 


3655 


7.172 


5.414 


5.543 


5.391 


1.785(5] 


1.828(5] 


1.777(5] 


0.3866 


0.3861 


0.3862 


3 D a -)■ 3 P!j 


3260 


0.485 


0.365 


0.374 


0.364 


8.541(3] 


8.750(3] 


8.510(3] 


0.0186 


0.0185 


0.0185 


^2b<a Aab 












4.623[5] 


4.722(5] 


4.602(5] 








3 D 2 -> 3 Pf 


3714 






16.605 


16.149 




3.448(5] 


3.354(5] 






0.8058 


3 D 2 -> 3 P§ 


3320 






5.602 


5.449 




0.831(5] 


0.808(5] 






0.1942 


^2b<a ^ab 














4.279(5] 


4.162(5] 








3 D 3 -> 3 P° 2 


3421 






31.519 


30.655 




3.652(5] 


3.552(5] 









TABLE VII: The lifetimes of the 5s4d 3 Dj states in ns. The 
last three rows give term-averaged 3 D lifetime. 





MBPT 


All 


Recomm. 


r(5s4d 3 Di) 
r(5s4d 3 D 2 ) 
r(5s4d 3 D 3 ) 


2163 


2113 
2337 
2738 


2171(24) 
2403(27) 
2816(31) 


r(5s4d 3 D) 
Expt. [15] 
Expt. [14] 




2453 
2900(200) 
2500(200) 


2522(28) 


by 








Ft (y) p 


4tt 3 32tt 5 


32tt 7 512tt 9 


(4) 


" 45jy + 189?/ 


f 45y 5 ' 99jy 7 ' 



The Eq. Q can be expressed in terms of the dc polar- 
izability a g (u) = 0) of the state g as [Tl| 

AE g = -^-(anf(k B T)±a g (0) + A£f°. (5) 
15 

Here A£^ yn is determined as 

and r\ represents a "dynamic" fractional correction to the 
total shift that reflects the averaging of the frequency 
dependence of the polarizability over the frequency of 
the blackbody radiation spectrum. Corresponding shift 
in the clock transition frequency, Ai/fpo_i So — (A£"^" — 

AEt^)/h, is referred to as dynamic correction to the 
BBR shift. The quantity 77 can be approximated by [Til ] 

80 7T 2 

63 (2J g + 1) a g {0jk B T 
^ \(n\\D\\ g )f ( 21vr 2 336tt 4 \ 



Dynamic corrections to the BBR shift of the 5s 2 1 So — 
5s5p 3 Pg clock transition in Sr at T = 300 K (in 
Hz) are given in Table [Vj The dynamic correction to 
the BBR shift of the 3 Pq level is dominated by the 
contribution from the 5s5p 3 Pq — 5s4d 3 Di transition, 
which contributes 98.2% of the total. Our final result 
—0.1492(16) Hz is in excellent agreement with recent 
value -0.1477(23) Hz of Ref. [3. 

Our result enables us to propose an approach for fur- 
ther reduction of the uncertainty in the BBR shift: a 
measurement of the hsAd 3 Di lifetime with 0.5% uncer- 
tainty would provide the value of the BBR shift in Sr 
clock that is accurate to about 0.5%, which would be a 
factor of 2 improvement in the uncertainty stated here. 
Such a determination assumes accurate knowledge of the 
branching ratios. 

The 5s4d 3 Di level decays to all three 5s5p 3 Pq.i 2 
states, but the branching ratio to the 3 P2 level is very 
small. The lifetime of a state a is calculated as 

1 

^2b<a Aab 

The El transition rates A a b are calculated using 
_ 2.02613 x 10 18 S _j 

A ab - ^3 2JaTT 8 ' 

where A is the wavelength of the transition in A and S is 
the line strength. 

We find that the branching ratios are essentially inde- 
pendent of the correlation corrections to the matrix ele- 
ments. We note that line strength ratios are close to the 
non-relativistic ones (5/9, 5/12, 1/36), with the differ- 
ences being -0.23%, +0.22%, and 1.4% for the 3 Di- 3 Pj 
transitions, respectively. We illustrate this point in Ta- 
ble [VH where we list the relevant energies, line strengths 
S 1 , transition rates A, and branching ratios in the CI, 
CI+MBPT, and CI+all-order approximations. 
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The RPA corrections to the matrix elements are in- 
cluded in all cases. We used experimental energies in 
all calculations for consistency. We find that the dif- 
ference in the CI, CI+MBPT, and Cl+all-order branch- 
ing ratio results is less than 0.1%. Since all the other 
corrections are small, their uncertainties should be even 
smaller. As a result, the accuracy of our branching ra- 
tios should be better than 0.2%. The recommended val- 
ues for the 5s5p 3 P° 2 — 5sAd 3 Di matrix elements are 
obtained using the recommended matrix element for the 
5s5p 3 Pq — 5s4d 3 D! transition and Cl+all-order branch- 
ing ratios. The recommended values for the transition 
rates and the hsAd 3 Di lifetime, 2172(24) ns, are ob- 
tained using the recommended values of the matrix ele- 
ments and experimental energies. We also list the recom- 
mended values 3 D2, 3 D3, and term-averaged 3 D lifetimes 
in Table PVTO The 3 D term-averaged lifetime is compared 
with experiment [bH . fisj ] . 

VII. CONCLUSION 

We have evaluated the static and dynamic polarizabil- 
ities of the 5s 2 1 So and 5s5p 3 Pg states of Sr and ex- 
plained the discrepancy between the recent experimental 



5s 2 1 So — 5s5p 3 Pg dc Stark shift measurement [l3[ and 
the earlier theoretical result Our theoretical value 

for the dc Stark shift of the clock transition, 247.5 a.u., is 
in excellent agreement with the experimental result. We 
have provided the recommended values of the matrix el- 
ements for transitions that give dominant contributions 
to the clock Stark shift and evaluated their uncertain- 
ties. We evaluated the dynamic correction to the BBR 
shift of the % - 3 Pg clock transition at 300 K to be 
-0.1492(16) Hz and proposed an approach for further re- 
duction of the uncertainty in the BBR shift. 
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